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The entanglement dynamics of an anisotropic two-qubit Heisenberg XYZ system in the presence 
of intrinsic decoherence is studied. The usefulness of such system for performance of the quantum 
teleportation protocol 7o and entanglement teleportation protocol Ti is also investigated. The results 
depend on the initial conditions and the parameters of the system. For the product and maximally 
entangled initial states, increasing the size of spin-orbit interaction parameter D amplifies the effects 
of dephasing and hence decreases the asymptotic entanglement and fidelity of teleportation. We 
show that the XY and XYZ Heisenberg systems provide a minimal resource entanglement, required 
for realizing efficient teleportation. Also, we find that for the some special cases there are some 
maximally entangled states which are immune to intrinsic decoherence. Therefore, it is possible to 
perform the quantum teleportation protocol 7o and the entanglement teleportation 71 with perfect 
5^ ' quality by choosing a proper set of parameters and employing one of these maximally entangled 

robust states as initial state of the resource. 
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I. INTRODUCTION 

Entanglement is a central theme in quantum information processing as a uniquely quantum mechanical resource 
that plays a key role in many of the most interesting applications of quantum computation and quantum information 
PHI]. Thus a great deal of efforts have been devoted to study and characterize entanglement in the recent years . 
The central task of quantum information theory is to characterize and quantify entanglement of a given system. A 
pure state of pair of quantum systems is called entangled if it is unfactorizable, e.g. singlet state of two half-spin 
7— I ' system. A mixed state p of a bipartite system is said to be separable or classically correlated if it can be expressed as 
a convex combination of uncorrelated states pa and pb of each subsystems i.e. p = u>ip A ® p B such that uii > 
and u>i — 1, otherwise p is entangled @, |j|- Many measures of entanglement have been introduced and analyzed 
[l|-[3j]. Here we use the negativity as a measure of entanglement. For the C 2 <g> C 2 bipartite systems, the negativity of 
C*~> a state p is defined as 

o ; 

E N (p) = 2max{-A min ,0}, (1) 

where \ m i n is the minimum eigenvalue of p TA , and Ta denotes the partial transpose with respect to the part A of the 
bipartite system. 

It is well known that a two-qubit entangled system can be used to perform the quantum teleportation protocols 
0. The pioneering authors of quantum infirmation theory have showed that the mixed quantum channels which 
allow to transfer the quantum information with fidelity larger than | are worthwhile 6]. By using the isomorphism 
between quantum channels and a class of bipartite states and twirling operations, Horodecki et al. have shown that 



the optimal fidelity of teleportation for a bipartite state acting on C 2 ® C 2 [7j is & mo .x = djr ^\ rl where J 7 , 



is 



the fully entangled fraction of the resource. Then, Bowen and Bose have shown that in the standard teleportation 
protocol Tq with an arbitrary mixed state resource the teleportation process can be considered as a general depolarizing 
channel with the probabilities given by the maximally entangled component of the resource and without additional 
twirling operations Using the property of the linearity of teleportation process, Lee and Kim Q have shown that 
quantum teleportation preserves the nature of quantum correlation in the unknown entangled states if the channel 
is quantum mechanically correlated. They have considered entanglement teleportation of an entangled state via 
two independent, equally entangled, noisy channel, represented by Werner states. In this teleportation protocol 71, 
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the joint measurement is decomposable into two independent Bell measurements and the unitary operation is also 
decomposable into two local Pauli rotations. In other word, T\ is a straightforward generalization of the standard 
teleportation protocol 7o just doubling the setup. Lee and Kim found that the quantum correlation of the input state 
is lost during the teleportation even the channel has nonzero entanglement. They also found that in order to achieve 
a favorite teleportation fidelity, the quantum channel should possess a minimal entanglement. Hence, in comparison 
with quantum teleportation, entanglement teleportation imposes more stringent conditions on the quantum channel 

Unfortunately, decoherence destroys the quantumness of the system and hence will decrease the useful entanglement 
between the parts of the system [ToMl2i | . There are several approaches to consider the decoherence and solve the 
quantum to classic transition problem. One of these approaches is based on modifying the Schrodinger equation 
in such a way that the quantum coherence is automatically destroyed as the system evolves. This mechanism is 
called "intrinsic decoherence" and has been studied in the framework of several models (see [l3[ and references 
therein). In particular, Milburn has proposed a simple modification of the standard quantum mechanics based on the 
assumption that for sufficiently short time steps the system evolution is governed by a stochastic sequence of identical 
unitary transformation rather than continuous unitary evolution [bH ] . This assumption leads to a modification of the 
Schrodinger equation which includes a term corresponding to the decay of quantum coherence in the energy basis. 
Using a "Poisson model" for stochastic time step, Milburn obtained the following dynamical master equation in the 
first order approximation 

±p(t) = -i[H,p]-±-[H,[H,p(t)]}, (2) 

where H is the Hamiltonian of the system, p(t) indicates the state of the system and 7 is the mean frequency of the 
unitary step and determines the rate of decoherence (T3 |. In the limit 7 — > 00 the Schrodinger's equation is recovered. 
Note that, in this mechanism of decoherence, the decay of quantum coherence is a result of phase relaxation process, 
so in the following we will only deal with dephasing processes without the usual energy dissipation associated with 
normal decay. The first order correction to the equation ([2]) leads to diagonalization of the density operator in the 
energy eigenstate basis, 

^ t (e'\p(t)\e) = ^(e - e'){e'\p{t)\e) - i-( £ - e'f '(s'\p(t)\e). (3) 

Note that the rate of diagonalization (dephasing) in the energy basis depends on the square of the energy separation 
of the superposed states. Thus the coherence between states that are widely separated in energy, decays rapidly. A 
formal solution of the Milburn's dynamical master equation @ can be expressed as [l3[ 



p(t) = ^2M k (t)p(0)Ml(t), (4) 

where 



fc=0 



M fe (t) = J ^H k exp(-iHt) exp(-^-H 2 ). (5) 
V re! 27 

It is evident from the master equation ^ that the state of the system remains constant in time, if the initial state of 
the system commutes with H. Thus all density matrices which can be written as a classical mixture of the eigenstates 
\ipi) of the Hamiltonian, i.e. p = ^2 Pi\ipi) (ipi\ with ^Pi — 1, are immune to intrinsic decoherence. One important 

i i 

case is the thermal state of the system with pi = t ^ e -/u^ , where T = is the temperature and ks is the Boltzman 
constant [l5j |. In other words, the set of all thermal states with different temperatures span an intrinsic decoherence 
free subspace. A decoherence free subspace is a Hilbert space such that each state of this space is immune to 
decoherence (lo| . 

The effects of intrinsic decoherence on the dynamics of entanglement and qua ntum teleportation and entanglement 
teleportation of Heisenberg systems have been studied in a number of works [16l-[l9j . For example Ye Yeo in Ref. [l6[ 
has shown that in an anisotropic two-qubit XY Heisenberg system, the nonzero thermal entanglement produced by 
adjusting the external magnetic field beyond some critical strength is a useful resource for teleportation via 7o and 7i 



protocols. The authors of Ref. [17[ have shown that, adjusting the magnetic field can reduce the effects of the intrinsic 
decoherence and accordingly one can obtains the ideal fidelity of teleportation via XYZ Heisenberg systems. Also, the 
results of Ref. [18j showed that for an initial pure state of the resource, which is the entangled state of a two-qubit 
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XXZ Heisenberg chain, an inhomogeneous magnetic field can reduce the effects of intrinsic decoherence. Then the 
authors of Ref. [19| argued that if the initial state is an unentangled state, we can improve the fidelity of teleportation 
protocol 7o via two-qubit Heisenberg XXX systems in the absence of the magnetic field by introducing the spin-orbit 
(SO) interaction , arising from Dzyaloshinski-Moriya (DM) interaction. However, the dynamics of entanglement and 
entanglement teleportation of more involved spin systems has not been discussed, yet. 

In this paper, we investigate the influence of the intrinsic decoherence (dephasing) on the entanglement dynamics 
and teleportation scheme of a two qubit anisotropic Heisenberg XYZ system under the influence of an inhomogeneous 
magnetic field and in the presence of SO interaction. This system is suitable for modelling of a system which is 
realized by the spin of two electrons confined in two coupled quantum dots [20l - [2^ |. Because of weak vertical or 
lateral confinement, electrons can tunnel from one dot to the other and spin-spin and spin-orbit interactions between 
the two qubits exist. In summary, we show that the dynamical and asymptotical behavior of the entanglement, the 
quality of the quantum teleportation and the entanglement teleportation and also the entanglement of the replica 
state, are dependent on the initial conditions and the parameters of the model. We discuss the problem for some 
special initial states and investigate the role of the parameters of model (such as coupling coefficients J M , magnetic 
field B, inhomogeneity of magnetic field b, SO interaction parameter D,...) on the entanglement properties of the 
system. The results show that for the product and maximally entangled initial states, the asymptotic value of the 
entanglement decreases with D. The fidelity of the teleportation approaches | form above for large values of D, for 
both product and maximally entangled initial states of the resource. Furthermore the results show that, it is possible 
to perform the teleportation protocols To and 71 with perfect quality in the XY and XYZ Heisenberg systems if the 
initial conditions and system parameters are set properly. 

The paper is organized as follows. In Sec. II, we introduce the Hamiltonian of a Heisenberg system under the 
influence of an inhomogeneous magnetic field with taking into account the SO interaction. For given initial states, 
the density matrix of the system at a later time is derived exactly by solving the Milburn's dynamical equation. The 
effects of initial conditions and system parameters on the dynamics of entanglement, as measured by negativity, is also 
studied in this section. The quantum teleportation and entanglement teleportation processes via the above system 
are investigated in the subsections II-A and II-B. Finally, in Sec. Ill a discussion concludes the paper. 



II. THEORETICAL TREATMENT 



The Hamiltonian of a two-qubit anisotropic Heisenberg XYZ-model in the presence of inhomogeneous magnetic 
field and spin-orbit interaction is defined by (l5j 

H = \{J X o\ol + J y a\a\ + J z o 7 {a z 2 + B x • <r x + B 2 • <x 2 

+ D • (<7i X <T 2 ) + 6 «Ti • T • <T 2 ), (6) 

where <Tj = (rrj, crjf, crj) is the vector of Pauli matrices, Bj (j = 1, 2) is the magnetic field on site j, J M (/i = x, y, z) are 
the real coupling coefficients (the chain is anti-ferromagnetic (AFM) for J M > and ferromagnetic (FM) for J M < 0) 
and D is Dzyaloshinski-Moriya vector, which is of first order in spin-orbit coupling and is proportional to the coupling 
coefficients (J M ) and T is a symmetric tensor which is of second order in spin-orbit coupling [23l - [26| . For simplicity, 
we assume Bj = Bj z such that B\ = B + b and B 2 = B — b, where b indicates the amount of inhomogeneity of 
magnetic field. If we take D = J Z D z and ignore the second order spin-orbit coupling, then the above Hamiltonian 
can be written as [30l |: 

H = JxiPx G t +a-^o-^) + (J + iJ x D)af<T2 + (J - iJ z D)a^a^ 

+ y^of + (— K + (— g— K. (7) 



where J := v , is the mean coupling coefficient in the XY-plane, x : — 'x.+ J j y specifies the amount of anisotropy 

in the XY-plane (partial anisotropy, — 1 < % < 1) and a = h( 
spectrum of H is easily obtained as 



in the XY-plane (partial anisotropy, — 1 < % < 1) and o ± — \{a x ± icr v ) are lowering and raising operators. The 



F|V ± )=£i,2|V ± ), 

(8) 

/J|E±) = £ 3 ,4|S ± ) , 
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where the eigenstates expressed in the standard basis {|00), |01), |10), |11)} are 

1^)= ^((7^)101) +|10», 



|E±)=M ± ((^2)|00> + |11}) 



(9) 

Mnn\ xln\) 

Jx 

with the eigenvalues 

(10) 

£3,4 = ± ) 



respectively. In the above equations N ± = 1 and M — , 1 are the normalization constants. 

fl | ("±Q 2 ~ / 1+ (S±2)2 
V 1+ ^ + (./,D)^ V JX 

Here we have defined, £ := ^/6 2 + J 2 + (J Z D) 2 and 77 := ^ B 2 + (Jx) 2 , for later convenience. 

According to Eq.(|4|) it is easy to show that, the time evolution of the density operator p(t) for the above mentioned 
two qubit Heisenberg system which is initially in the state p(0), under intrinsic decoherence is given by 

p(t)= ^2 exp[ — -(e m - e n ) 2 - i(e m - e n )t] (<f> m \p(0)\^> n ) \<j> m ){<j>n\, (11) 

m,n— 1 

where eigenenergies e m . n and the corresponding eigenstates 101,2} = and 103,4) = l^*) are given in Eqs.© and 
(jT0|) and 7 is the phase decoherence rate. 

In the following we will examine the evolution of entanglement under intrinsic decoherence of a class of bipartite 
density matrices having the standard form 

P (o) = ^ + |oo)(oo| + M-|n)(ii| + K|oo)(ii| + |ii)(oo|) 

+ wi|01)(01| + ™ 2 |10)(10| + z|01)(10| + z*|10)(01|, (12) 

which is called X states class and arises naturally in a wide variety of physical situations. If the initial state belongs to 
the set of X states (fT2|) . then Eq. |(5J) guarantees that p(t) given by Eq. ([TTj) also belongs to the same set. Therefore, 
the only non- vanishing components of the density matrix in the standard basis are 

Pll (t) = ^{2B 2 +(j X ) 2 (i+$(t))] + ^[(Jx) 2 (i-m)}-^{BJx(i-m)}, 



2rj 2 2i] 2 t] 2 

/-( / ) = ||[2& 2 + (J 2 + (J^) 2 )(l + <*>'(*))] + ||[(J 2 + (J z D f)(l <!>'(/ ) I 

[JL-[(J-iJ z D)b(l-*'(t))} + C.C], 



P33(t) = ^[(J 2 + {J z Df){l^^\t))] + ^[2b 2 + {J 2 + (J z Df)(l + ^{t))] 



2£ 2 

[^[(J-iJ z D)b(l -*'(*))] + C.C\, 
P 44 (<) = 0[(J X ) 2 (1 -$(*))] + ^[2S 2 + (J X ) 2 (l + $(t))] + ^[BJx(l-$(t))], 



2rj 2 2?y 2 rf 

Arj* 2rj z r/ z 



P2 3 (t) = ^[(J + iJ z D)b(-l + *'(t))} + ^[(J + iJ z D)b(l-*'(t))} 

+ ^2 [( J 2 + ( J ^) 2 ) + e Wl + |*K J + iJ * D ?{i *'(*)))] » 

P4i(*) = , 

032 (*) = P23(t)*. (13) 



Here, we have defined $(t) := e^ 2 ^ 7 * cos2r/i, *(t) := [cos2r?t - § sin2?7t]e~ 2 '' 2 '»'*, $'(*) := e" 2 ^ 7 * cos 2£f, := 

b 



[cos2£i - f sin2£i]e _2? 7 * and 6(t) := [(£ 2 + & 2 )cos2£t - 2^£sin2^]e" 2 « 7 *. For asymptotically large times, all 
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FIG. 1: (Color online) Dynamical and asymptotical behavior of the negativity vs. the model parameters and for the product 
initial states, |00) (left graph) and |01) (right graph). The parameters of the model are chosen to be J = 1, \ — 0.9, J z = 0.5 
and 6 = 1. We have set D = and 7 = 0.09 for the left and B = 3 and 7 = 0.02 for the right graph. All parameters are 
dimensionless. 



of these functions vanish, and hence the state of the system at asymptotically large time limit p° 
easily. Knowing the density matrix pit), we can calculate the entanglement by negativity: 



where 



Nip) = max{-2min{Ai, A 2 , A 3 , A 4 },0} 



can be obtained 



(14) 



Ai, 2 = -((wi + w 2 ) ± y/((>22(t) - pasit)) 2 + 4|puW| 2 ) 



A 3 ,4 - -(( M + + M-) ± VOM*) - Mi)) 2 + 4|p 23 (t)| 2 , 



(15) 



are eigenvalues of the partially transposed matrix p TA it). The negativity is a function of the model parameters and 
the initial conditions. Figs. [1] and [5] depict the time variation and asymptotical behavior of the negativity for product 
and maximally entangled initial states, respectively. The results show that, the entanglement reaches a steady state 
value after some coherent oscillations for the times greater than -. The size of this steady state value, 7V°°, depends 
on the initial conditions and the parameters of the model. Fig. [T]shows that for the product initial state |?/>(0)) = 1 00) , 
the function N°° increases with magnetic field for B < 1 and then decreases for B > 1. Hence, for the case of B < 1 
and IV'(O)} = |00), we can suppress the effects of dephasing by increasing B. This figure also shows that, for the case 
of product initial state | -0 (0) } = |01), iV°° is a decreasing function of D, for all values of D. Fig. [5] shows that, for 
the case of maximally entangled initial states |^(0)) = ^° ^ 11 ^ and |^(0)) = ^j^ 1 ; the function TV 00 decreases 
with B and D, respectively. This figure also shows that, in the case of B — 0, the state |^(0)) = hag robust 

entanglement. Also the state |"0(O)) = ^ 01 ^ 10 ^ has robust entanglement, when D = b = 0. 
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FIG. 2: (Color online) Dynamical and asymptotical behavior of the negativity vs. the model parameters and for the maximally 
entangled initial states, l 00 '^ 11 ' (left graph) and ^ 01 ^ 10 ^ (right graph). The parameters are the same as in Fig. [T] 



A. Quantum Teleportation 



According to the results of Bowen and Bose the standard teleportation protocol 7o, when used with two-qubit 
mixed state of the Heisenberg spin chain p(t) as a resource, acts as a generalized depolarizing channel A^[p in }. In 
the standard teleportation protocol an input state is destroyed and its replica (output) state appears at remote place 
after applying a local measurement and unitary transformation in the form of linear operators. We consider as input 
an arbitrary pure state \tpi n ) = cos ||0) + sin ||1)(0 < < n, < (f> < 2n). The output (replica) state, p out , can 



be obtained by applying joint measurement and local unitary transformation on the input state p^ 
state is given by 



Pout 



4 

£ 



Pfi cr P Pin O^, 



Thus the output 



(16) 



where p — 0, x, y, z (cr° — I), p u — Tr\E^ L p(tj\ represents the probabilities given by the maximally entangled fraction of 
the resource p(t). Here E° = |*-)(*-|, E 1 = |$")($-|, E 2 = |$+)(<I> + and E 3 = |*+)(*+| where = (|01> ^ ln>) 

and = (|00) ^ n)) are the Bell states. 

The quality of the teleportation is characterized by the concept of fidelity. The maximal teleportation fidelity 
achievable in the standard teleportation protocol 7o is given by Q 



.(A#>) 



2T{t) + 1 



(17) 



where JF{t) = max {pA is the fully entangled fraction of the resource. For our model, the probabilities, p„s can 

u— 0,1,2,3 

be written as 



Po,3 = +w 2 ) ± Ue[p 2 3{t)], 

Pi, 2 = \{p+ + P-) ± ^e[p 14 (t)}. 



(18) 
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Therefore, the maximum fidelity depends on both the initial conditions of the quantum channel and the parameters 
of the model. In the following let's examine some important entangled and product initial states for the channel: 
i) 1^(0)) channel = -/g i m this case we have 101,2 = 5, z = ±5 and (i±=v = and hence 



wAg>H4 J2+(t2+ r m( " . 



For the asymptotically large times vanishes and we have 

This equation states that, the maximum fidelity achievable at large time limit is always greater than | i.e. this 
channel is superior to the classical channels. The function <&^ ax is minimum (Q^ax = §), if the interaction on the 
resource is Ising type in the z direction (i.e. J=0) and reaches its maximum {<&%° lax — 1) f° r the case of D = 6 = 0. 
Note that in the later case the state of the channel is a maximally entangled state (see Fig. [5]). 

ii) \ip(0)} channel = ^^3^' i- e - M± = \ i v = ±^ an d ^1,2 = z = 0. The maximum fidelity achievable for this 
quantum channel is 



and hence, 

*™«* = *ma*(AO = i + ^) 2 ' ( 22 ) 



I 4- hi** 

10 3 3 77 



which means that, the XY and XYZ chains (J, \ ^ 0) are more useful resources for performance of the teleportation 
protocol To- In this case, we have Q^ax = 1 for £? = 0, this result is compatible with the results of Fig. J2|). 
Hi) \ip(0)) channel = 1 01) (or 1 10} ) , in this case just Wx = 1 (^2 = 1) is nonzero and hence 



At the large time limit we can write 



*~ax = *r n ox(AC) = | + ^f- (24) 



and in this case, we can adjust the quality of quantum teleportation by changing b, J z and J. The asymptotic fidelity, 
^maxi tends to I from above for large values of b, thus in the case of J ^ our channel is superior to the classical 
communication. In this case, increasing \ Jz\, decreases the quality of teleportation, and hence the XY chain is more 
suitable than XYZ chain. There is no way to reach the value $^ a2 . = 1, since the equation \bJ\ = £ 2 has no real 
solution. 

iv) \ip(0)) channel = |00) (or |11)), i.e. /_t + = 1 (/i_ = 1). In this case we have, 

w^)-l + j|j y m)| . <*> 

and for asymptotically large time we have 

$~a* = ^a,(AC) = \ + (26) 

According to the Eq. (|26l) . the presence of anisotropy in XY-plane (x 7^ 0) provides the desirable fidelity ($°° > |). 
The fidelity cannot take the maximum value $^ aa . = 1, because the equation -B|Jx| = ?/ 2 has no solution in the 
domain of real numbers. 
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FIG. 3: (Color online) Dynamical and asymptotical behavior of the output entanglement of the entanglement teleportation 
protocol 71 vs. the model parameters and for the product initial states, \tp(Q)) channel = 1 00) (left graph) and \ip(0)) channel = 1 01) 
(right graph). The input state is considered a maximally entangled state Ni n = 1. The parameters are the same as in Fig. Q] 



B. Entanglement Teleportation 

In this section, we consider Lee and Kim's teleportation protocol 71 and use two copies of the above two-qubit 
state, p(t) (g) p(t), as resource [9]. We consider the pure state \ip in ) = cos § 1 10) + sin02|Ol) (0 < 9 < vr, < <j> < 2tt) 
as the input state. The negativity associated with the input state, pi n = |V'm)(V'm| is ^(pin) — Ni n = sm9. By 
generalizing Eq. (|16[) the replica (output) state p ou t can be written as 



Pout 



k p ^® p{t) Pm = J> pi > M ® T^An^ ® <r v ), (27) 



where p = 0,x,y,z, p av = p a Pv By considering the two-qubit spin system as a quantum channel, the state of the 
channel is given by the equation (1131) and hence one can obtain p out as 

Pout = a(|00)(00| + |ll)(ll|) + K (<)(|00)(ll| + |ll)(00|) 

+ a '|01)(01| + c'(t)|01)(10|+c'*(<)|10)(01|+6'|10)(10| (28) 

where 



a 


= (lOl + w 2 ){p + + p ), 




«(f) 


= A^e[p 23 {t)}^e[ Pli {t)} cos 




a' 


= (p+ + p-) 2 cos 2 | + 


-w 2 ) 2 sin 2 f, 


b' 


= (wi + w 2 ) 2 cos 2 § + (/i+ 4 


~p-) 2 sm 2 ^, 


c'(t) 


= 2 e -^((5Re[p 23 W]) 2 + e 2 ^ 


{fte{p 14 (t)}) 2 ) sin 



(29) 

Now, we can determine the negativity of the output state as 

N out = N(p out ) = max{-2min{A' 1 , A 2 , A 3 , A^},0} (30) 
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FIG. 4: (Color online) Dynamical and asymptotical behavior of the output entanglement of the entanglement teleportation 
protocol 71 vs. the model parameters and for the maximally entangled initial states, \ip(0)) channel = 00 ^ 1:L ^ (left graph) and 
IV'(O)) channel = 01 y^ 10 ^ (right graph). The input state is considered a maximally entangled state Ni„ = 1. The parameters 
are the same as in Fig. [T] 



where, 

V 1)2 - \{{a' + tf)±y/{p! -V)*+4\n{t)\*) 

\' 3A = a±\c'(t)\, (31) 

are the eigenvalues of p^£t (t) . The function N out is dependent on the entanglement of the input state N in and the 
entanglement of the resource N c h a nnei (which is determined by the initial condition and the parameters of the channel) . 
The dynamical and asymptotical behavior of the negativity of output state of the entanglement teleportation protocol 
7i are illustrated in Figs. [3] and 0] for the product and maximally entangled initial states of the resource, respectively. 
In these figures, we assume that the input state is a maximally entangled state Ni n = 1. The results are compatible 
with the results of Figs. [T] and [5] By analyzing these figures one can find that, more entangled channel state is, 
more input entanglement is preserved. Figure [5] depicts the behavior of the N™ t versus the spin-orbit parameter D 
(which indicates the entanglement of the channel) and iVj„, for a given set of other parameters. Increasing the value 
of D causes N^ annel to decrease (see Figs. Q]aJid|2} and hence, as the Fig. [3] shows, the entanglement of asymptotic 
output state, decreases as D increases, for a fixed value of iVj n . Fig. [5] also shows that, for a fixed value of D (i.e. 
fixed N^ anne[ ), N™ t is an increasing function of N in . This means that as the input entanglement increases, a more 
entangled quantum channel is required to realize efficient entanglement teleportation. 
The fidelity between pi n and p out in terms of input negativity (JVj n ) is obtained as (27j 

F(p in , Pou f,t) = \(lpin\poutt\^in)\=fl(t) +f 2 (t) Nf n , (32) 

where f t (t) = (wi+m 2 ) 2 and f 2 (t) = \ - (w x + w 2 ) + 2((Ke[pi 4 (i)]) 2 cos 2<f>+ {^e[p2z{t)]) 2 ) . For <j> = 0, these functions 
depend only on the parameters of the channel. This formula has been also reported in Ref. [|| by Kim and Lee, but 
contrary to the Werner states, / 2 (i) can be a positive number for Heisenberg chains. This means that, there exists a 
channel which teleports more entangled initial states with more fidelity, but it should be noted that, if we choose the 
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FIG. 5: ((Color online) Asymptotic negativity of the output state of the entanglement teleportation protocol Ti vs. D and the 
entanglement of the input state, Ni n . The parameters are the same as in Fig. [1] 



parameters of the channel such that f 2 (t) > then f x (t) decreases and ultimately F(pi n , pout', t) becomes smaller than 
g, which means that the entanglement teleportation of mixed states is inferior to classical communication. Thus, to 
obtain the same proper fidelity, more entangled channels are needed for more entangled initial states. 

The average fidelity Fa is another useful concept for characterizing the quality of teleportation. The average fidelity 
Fa of teleportation can be obtained by averaging F(pi n , pout) t) over all possible initial states 



F A {t) 



J 27r d(j> J* F(p in , p out ; t) sin OdO 
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±(2( Wl + w 2 f + (p + + ,i_) 2 + 4(5J e [ P23 (i)]) 2 ) . (33) 



The function Fa(£) depends on the initial conditions and parameters of the channel. According to this formula the 
asymptotic value of the average fidelity of entanglement teleportation becomes = |(1 + ^(j) 4 ), if the channel is 

initially in the state, \ip(0)) channel = ■ Fig. |6]gives a plot of F^ in terms of the parameters D and b (which 

determine N c hannel), in this case. The figure shows that, for fixed values of the other parameters and b, F^ decreases 
as D increases (or equally, N c h a nnei decreases), such that for the large values of D, F™ approaches | from above. We 
can achieve perfect entanglement teleportation (FT = 1) in the case of D = b = 0. For the case of product initial 
state, \tp(0)) channel = |01) of the resource, we have F™ = |(1 + ^(|r) 2 ) which tends to | from above for large values 
of D, too. But never reaches the value 1 for any choice of the parameters. 



III. DISCUSSION 



The effects of dephasing due to intrinsic decoherence on the entanglement dynamics of an anisotropic two-qubit 
Heisenberg XYZ system in the presence of an inhomogeneous magnetic field and SO interaction, are investigated. The 
usefulness of such systems for performance of the quantum teleportation and entanglement teleportation protocols 
are also studied. Intrinsic decoherence destroys the quantum coherence (and hence quantum entanglement) of the 
system as the system evolves. For the case of noninteracting qubits dephasing processes kill the quantum correlations 
(entanglement) of the system at a finite time and hence entanglement sudden death (ESD) phenomenon occurs (i.e. 
entanglement vanishes faster than local coherence of the system [28|, The results of this paper shows that, for 

interacting qubits, dephasing induced by intrinsic decoherence is competing with inter-qubit interaction terms to 
create a steady state level of entanglement after some coherence oscillation and hence the entanglement of the system 
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FIG. 6: (Color online) Asymptotic mean fidelity of the entanglement teleportation protocol 71 vs. D and b for maximally 
initial state, l 01 ^! 10 ^ of the resource. The parameters are the same as in Fig. Q] and the value of b = is chosen for the right 
graph. 



reaches a stationary value, asymptotically. The dynamical and asymptotical behavior of the entanglement depends 
on the initial conditions and the system parameters. Indeed, the effects of dephasing can be amplified or weakened 
by adjusting the parameters of the model and initial conditions. We show that for the product and maximally 
entangled initial states, the asymptotic value of the entanglement decreases as D increases. This is because due to 
hermiticity of the Hamiltonian, we can express the state of the system as a superposition of the energy eigenstates. 
Increasing D, increases the energy separation of the superposed states (see Eq. (fT0|)) and hence amplifies the effects 
of dephasing (see Eq. Consequently, for both product and maximally entangled initial states of the resource, 

the fidelity of teleportation approaches | form above for large values of D, this is the maximum fidelity for classical 
communication of a quantum state. Furthermore, our results show that, for product initial states and a specific 
interval of the magnetic field B, the asymptotic entanglement (and hence the fidelity of the teleportation) can be 
enhanced by increasing B. We also have argued that, a minimal entanglement of the resource is required to realize 
efficient entanglement teleportation. The results also show that, the XY and XYZ Heiscnbcrg interaction can provide 
this minimal entanglement for the channel state. 

We have also found that, the thermal state of the system, px = j^rjrpm is immune to intrinsic decoherence and 
hence it has robust entanglement with respect to the intrinsic decoherence. In the absence of magnetic field (B = 0), 
the maximally entangled initial state, | -0(0) ) = l 00 ^! 11 ^ , are immune to intrinsic decoherence and consequently, have 

robust entanglement. The same result is also true for the maximally entangled initial states \i{)(0)) = ^ 01 ^ 10 ^ in the 
absence of SO interaction (D = 0) and for homogeneous magnetic field (b — 0). Therefore choosing a proper set of 
parameters and employing one of these robust states as initial state of the resource, enable us to perform the quantum 
teleportation protocol 7o and the entanglement teleportation 71 with perfect quality (& m ax = Fa = 1)- 
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